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Abstract 

The a-cluster model is based on two assumptions that the proton-neutron pair interactions 
are responsible for the adherence between a-clusters and that NN-interactions in the a-cluster 
I structure are isospin independent. It allows one to estimate the Coulomb energy and the short 

. range inter-cluster bond energy in dependence on the number of clusters. Charge radii are 

' calculated on the number of a-clusters too. Unlike the Weizsacker formula in this model the 

. binding energies of alpha-clusters and excess neutrons are estimated separately. The calculated 

values are in a good agreement with the experimental data. 
. keywords: nuclear structure; alpha-cluster model; Coulomb energy; Surface tension energy, 

binding energy; charge radius 
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O ■ 1 Introduction 
m ; 

\^ \ The liquid drop model is explicitly formulated by two phenomenological formulas, the Weizsacker 

' formula [1] to calculate the nuclear binding energy and the formula constA^^^ to calculate the 

r~| . nuclear radii. It is known [2] that the idea of alpha clustering has developed together with the 

I , liquid drop model since the beginning of the nuclear physics. Nowadays the alpha-cluster model 

' is becoming more popular. Recently it has been shown that a large body of data on the a-decay 

2 ■ is described by a simple cluster model with using a notion of an a-core potential [3]. There are 

. . \ several attempts to incorporate the alpha-cluster model into the shell model [2,4]. To elucidate the 

I optical potential in nucleus-nucleus collisions an idea of bigger clusters like nucleon molecules of 

^ ' ^^C has been used [5]. Currently the idea of alpha-clustering is actively used in the microscopic 

studies of the nuclear structure of the light nuclei in the framework of the few-body approach [6] , 
as well as in the microscopic approaches to the effective nucleon- nucleon interaction [7] . Also there 
have been some studies recently to apply an alpha-cluster model to facilitate the description of 
some characteristics of nuclei like the a-particle separation energy [8] and the nuclear radii [9] for 
a wide range of nuclei. 

There is another strong evidence in favor of an alpha-cluster model of nuclear structure. It 
comes from the recent findings [10,11] that the radii of nuclei can be described independently of 
the number of excess neutrons by the formula R = const{Zy^^ fm. It can be considered as a 
consequence of a well known fact that the charge density in nuclei is almost constant (see, for 
example, Table 3.7 [12]). 

In this connection it is interesting to obtain the formulas to calculate the charge radii and 
binding energies of nuclei on the basis of the idea that the nuclear structure is defined by its 
elements, a-clusters. In such an approach, unlike the Weizsacker formula, the formula to calculate 
the nuclear binding energy should consist of two separate parts. One is responsible for the binding 
energy of a-clusters and the other should estimate the binding energy of excess neutrons. 

The core starts growing with Z from the nucleus with Z > W [10] together with the surface 
tension energy. The distance between the nearby clusters is determined by their charge radii, which 
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axe bigger than their matter radii. It was shown [10,11] that the /3-stable nuclei have a particular 
number of excess neutrons needed to fill in the difference in the volumes occupied by the charge and 
the matter of the nucleons of the a-clusters in the core. Then the radius of a /9-stablc isotope 
can be estimated by the volume occupied by the matter of the core and the volume of the charge 
of a few peripheral clusters. It has been shown that the condition Rm = R determines the narrow 
/3-stability path and its width. The binding energy of these excess neutrons can be calculated as a 
sum of the energies of nn-pairs. 

The hypothesis that the nra-pairs are placed in the core explains its high mass density. In 
addition, it provides an explanation [10] why the experimental values of the most abundant isotopes 
Rgxp well described by both functions, ~ A^/^ in the liquid drop model and ~ Z^^^ in the a- 
clustcr model. Then it is expected that the NN-interaction in the regular ct-cluster structure is 
different from that in excess neutron matter. 

It is approved by a recent investigation of the optical potential for the elastic channels in reaction 
^^Ca(p,n)^^Sc [13]. It has been found that the projectile proton is involved in the NN-interactions 
in p Ca channel somewhat different from NN-intcractions of the outgoing excess neutron in 
n Sc channel. The proton belongs to the a-cluster structure and the neutron is from the excess 
neutrons. To elucidate the isospin independence of NN-interactions between the nucleons of the 
a-cluster structure the same investigation should be done for the different reaction ^^Ca(p, n)^^Sc 
where the both proton and neutron belong to one pn-pair. 

The energy of the short-range inter-cluster bonds e^a is taken to be equal e^a = ^aa^~^aa where 
^aa denotes the nuclear force energy and e^^ denotes the Coulomb repulsion energy between two 
nearby clusters. The old term 'nuclear force energy' means here the energy of the attraction 
between nearby clusters due to the strong inter-nuclcon potential. The values of e^a = 2.425 
MeV, e^^^ = 4.350 MeV and = 1.925 McV have been found from analysis of the experimental 
data of the lightest nuclei with Z < 6 [14] under two assumptions that the proton-neutron (pn) 
pair interactions are responsible for the adherence between a-clusters and that the proton and the 
neutron belonging to one pair have equivalent single-particle nuclear bound state potentials, the 
EPN requirement [14,15]. The latter allows one to take the difference between the single-particle 
binding energies of the proton and the neutron of the last pair AEpn as the Coulomb energy of the 
proton's interaction with the other protons of the nucleus. 

It was shown [14,15] that the experimental binding energies .Egxp of the a-cluster nuclei with 
N = Z > 6, known for the nuclei with Z < 28, are well described by the formula E'^ = NaSa + 
3(-/V'a ^ 2)eaa where = Cffle = 28.296 MeV. The number of bonds 3{Na — 2) means that every 
new cluster brings three bonds with the nearby clusters. This idea was discussed before [2]. In the 
nuclei heavier than ^^O the long-range part of the Coulomb interaction of the last a-cluster with 
the remote clusters starts growing. It must be compensated by its surface tension energy E^. Then 
the empirical values of the surface tension energy of the last a-cluster E^ can be obtained from 
Eq. E^ = AEa — (e^ -|- Se^^,) where stands for the internal Coulomb energy of one a-cluster, 
which equals the Coulomb energy of the nucleus ^He = 0.764 MeV, and AE^ = AEpn of 
two pn-pairs belonging to one cluster [14,16,17]. 

The long-range part of the Coulomb energy of the last cluster's interaction with Na — 4 remote 
clusters is well approximated by expression 2{Z — 8)e^ /{1.2Rp) where Rp is the radius of the last 
proton position in the nucleus. Then it is expected that E^ = 2{Z — 8)e^ /{1.2Rp) = 71 -Rp, which 

determines Rp. To find the Coulomb radius Rc, the approximation E^ = JiRq = 72-^0^"^ is 
used. From the analysis of the empirical values of AE^ known for the nuclei with Z <22 with an 
accuracy of a few KeV the values of 71 and 72 have been obtained together with the formulas to 
calculate Rp and Rc on the number TV^. 

TliTis, the formula to calculate the binding energy with using the nuclear force energy of a- 
clustcrs = iVae^^ + 3(iV„ - 2)eS^^ the Coulomb energy E^ = 3/5{Zef/Rc and the surface 
tension energy 

^st ^ ^A^. ^st pi^g ^j^g binding energy of the excess neutrons -Eaa/" estimated as 
a sum of separation energies of nn-pairs describes the binding energy of the nuclei with > 5, 
with the errors close to those of the Weizsacker formula [16,17]. 
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It is evident that the model can be applied only for the nuclei with the number of excess neutrons 
AN > 0, so that every proton is coupled with a neutron. This set of nuclei includes all stable nuclei 
and the unstable nuclei around the /3-stability path. The model has been developed for the last few 
years and the main formulas have been published [10,11,14,15,16,17]. In the most recent papers 
[16,17] the surface tension energy was calculated with using the empirical values of for Z < 28, 

< 14. In this paper the value E^* is calculated by one formula for 9 < < 118. Besides, the 
phcnomenological formula to calculate the excess neutron binding energy has been extended for 
the cases of odd number of excess neutrons. 

2 Nuclear Binding Energy 

2.1 Nuclear force energy, Coulomb energy and surface tension energy 

The picture given in Fig. 1 shows the pn-pair bonds in the lightest nuclei. The distance between 
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Figure 1: A schematic illustration of two kinds of bonds between nearby pn-pairs, which define 
short-range inter-cluster bonds (solid line ellipse) and the bonds between the single pn-pair in the 
nuclei with odd Z and the pn-pairs of nearby a-clusters (dashed line ellipse). 

clusters is supposed to be constant for all nuclei except the nucleus ^Be. In this case one bond 
is not enough to keep two positively charged clusters close. Adding one more pn-pair causes two 
new bonds with the pairs of two alpha-clusters with the energy 2epnpn in the nucleus ■'^'^B, and the 
clusters get closer with the distance proper for the alpha-cluster liquid. The binding energy of the 
lightest nuclei E^ is calculated in accordance with Fig. 1 by the following equations: 

E\Ane) = £«, 

E (6Li) = S(y^ + Spnpn ~l~ ^pn-i 

E (lOB) = 2£'q; -|- £(x(x ~\~ '2iSpfipYi -\- Spfi-, 

E\12C) = 3£„ + 3£aa, 

where Saa denotes the binding energy between two nearby alpha-clusters, Spnpn denotes the binding 
energy between the single pn-pair and a pn-pair of the nearby a-clusters and £p„ stands for the 
binding energy of the proton and the ncTitron in the single pair. Solving these four equations 
with the values of E^ equal to the experimental values -Eexp [18] with an accuracy of 1 KeV (at 
experimental measurement error within 2-3 KeV) gives the values Sa = 28.296 MeV, £0,0 = 2.425 
MeV, epnpn = 2.037 MeV, Spn = 1.659 MeV. 

The experimental values of the other nuclei including ^^C are well approximated by means of 
the following equations for the nuclei with even Z with the mass number A: 

E^ = N^e^ + 3(iV„ - 2)£„„, (1) 

and for the nuclei with odd Zi = Z + 1 with the mass number Ai= A + 2 
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El — + 6£pnpn + £pn- (2) 

The energy of nucleus Ai is determined by the energy of the nucleus A plus the energy of 
six bonds between the single pn-pair and the six pn-pairs of the three nearby clusters plus the 
binding energy of the pair itself. The values of E^ are presented in Table 1 in comparison with the 
experimental values Eexp- 



Tabic 1. Binding energies calculated for the nuclei with N = Z. The values are given in MeV. 
Empirical values are given to three places of decimals and calculated ones are given to one place 

of decimals. 
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The other values in Table 1 have been defined in Introduction and are provided by the number 
of the corresponding formula. The values for even Z and for odd Zi nuclei, namely E^ and E^, 
E^ and , E"^""" and ^2r> ^''^and -Ef* are indicated as E^, E^"^", E^^ for the nuclei with Z with 
the mass number A. The 9th column titled as 'Ai?^(7)' gives the difference between E^xp and E'^ 
calculated by (7). The column with 'A£^^^(10)' contains the differences between Egxp and E^^ 
calculated by the Weizsacker formula(lO). The values indicated by '*' have been obtained by means 
of (17)/(18) and (21), see Section 3, with the accuracy of several tenth of MeV. 

The EPN requirement comes from isospin independence of NN-interactions in a-clustcrs. It 
provides a possibility to estimate the total Coulomb energy of a nucleus, E'-'' = X] ^Epn, where 
AEpn is the differences between the binding energies of the proton and neutron belonging to one 
pn-pair. If a sum of the values AEpn of two pn-pairs belonging to one a-cluster is presented as the 
Coulomb energy of the a-cluster AE^ = ^Epn, then 
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E^ = ^AEa + 5, (3) 
1 

where, the value of S comes from taking into account that two clusters in ^Be are at bigger distance 
than in ^°B, see Fig. 1. For the nucleus Ai the corresponding equation is 

= Y,I^E^ + 5 + AEpn. (4) 
1 

The values of AEp^ for the nuclei with N = Z are obtained by using the experimental binding 
energies taken from [18,19]. The binding energies of the nuclei with (Z — 1,N) and {Z,N — 1) 
needed for calculation of AEpn for the nuclei with 3 < Z < 8 and Z = 24, 25 have been calculated 
from the mass deficiency [19]. The energy of '^He 28.296 MeV and of ^He 27.338 MeV are taken 
from [18]. For the nuclei with Z = 23,26,27,28 the values of AEpn cannot be obtained by this 
procedure due to the lack of experimental data for the nuclei with {Z,N — 1). The values have 
been estimated in the framework of this model by equations (17)/(18) with using the values of the 
radius of the last proton position Rp/pi (21), see Section 3. 

The value of the Coulomb repulsion energy between two nearby clusters e^^ = 1.925 MeV has 
been found from the equation for the Coulomb energy of one cluster in AE^^C) = e^ + 2e'^^, 
where e'^ denotes the internal Coulomb energy of one cluster = Ef^^ = 0.764 MeV. This value 
perfectly fits the corresponding equation for ^^O, which is AEa(^^0) = + 36^^. The value of 
the Coulomb repulsion energy between the single pn-pair and one nearby alpha-cluster in case of 
odd Zi is found from the data for the nuclei ^Li and ^^N, e^^a = 1-001(6) MeV. 

One can estimate the Coulomb repulsion energy between two alpha clusters in the ^Bc as follows, 
£aa(^Be) = Epgg - 2eg = 1.887 MeV, which is less than e^^ by the value of 6 = 0.038 MeV. 

Eqs. (3) and (4) give the empirical values of the total Coulomb energy E'-^ obtained on the basis 
of the values known from the experimental data with an accuracy of a few KeV. Having obtained 
the total Coulomb energy of nuclei as a sum of AEpn of the pairs, one can easily test the validity 
of the hypothesis of the alpha-cluster structure on the nuclei with few clusters < 4, that is the 
nuclei where each cluster is in touch with everyone of the other clusters. In this case the Coulomb 
energy can be calculated by the number of clusters and their bonds. 

For ^Be according to Fig. 1. ^^(^Be) = 2e2+e^a{^Be). For and E^{^^C) = 3eg+3eg„ 
and E^(^^0) = 4eC + 6eg„. In case of ^Li E^'{^U) = e'^ + e^^^ and for ^^B E^'(^^B) = 2^^ + eC^ + 
AEpn. For ^^N £;'^(14N) = E'^{12C) + 3e^^^. For these nuclei the equations give values of E'^ in 
agreement with the values obtained by (3) and (4), see Table 1, with an accuracy of 1 KeV, for 
^Li with an accuracy of 6 KeV. It confirms the validity of the cluster structure of the nuclei with 
the quantities of e^Q, and e^^,, as well as the equations (3) and (4) are confirmed as the reliable 
formulas for obtaining the empirical values of the Coulomb energy for the nuclei with Z <22 with 
the sure accuracy of few tens KeV. It also approves the statement about isospin independence of 
NN-interaction. 

After adding the values of AEpn for the nuclei with Z = 26, 28 and Z = 29 estimated in the 
framework of the model with an accuracy of several tenth MeV (the accuracy is determined by the 
difference in the values obtained by the other possible way, by using (19)/(20), see section 3), the 
accuracy of the value of £''^ for ^^Cu is expected to be about 1-2 MeV. 

One can estimate the nuclear force energy of a-clusters from Eq. e^"'^ = + ~ 29.060 MeV 
and e'^l'' = e^a + = 2.425MeV + 1.925MeV = 4.350MeV. This consideration allows one to 
estimate eS^"(^Be) = (Eexp(^Be) - 2e„) + e'^^i^Be) = 1.792 MeV. 

The empirical values of nuclear force energy of a-clusters E^^"'^ for the lightest nuclei Z < 6, 
Na < 3, are obtained from E^xp = E'^^'^ + E'^ and they are presented in Table 1. For the nuclei 
with Z > 7, Na > 3, one gets 

rnnuc AT ^nuc i r>( kt o^^""c. ir^iuc jpnuc ■ «^nuc i _ /rx 

Ea - Na^a + ^V^a - ^J^aa > -^al - E + 'OSpnpn + %n, (p) 
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where Qe^npn + ^pn = ^^pnpn + 3ep„Q + ^pn = 16.884 MeV. 

The Coulomb energy for the nuclei with Na > 4 consists not only of the energies of bonds 
between nearby clusters. It must contain the long-range Coulomb interaction part. That fact that 
Eq. (l)/(2) has a good agreement with the experimental values allows one to suggest that the 
long-range part of the Coulomb energy E^^^" between the last cc-cluster and the Nq, — 4 remote 
clusters is compensated by the surface tension energy, which means that the E^^ ^ for the nuclei 
with Z >9. It gives an opportunity to obtain the empirical values E^^ from the values of AEa and 
AEpn, known with a good accuracy. 

The Coulomb energy of the last a-cluster for the even nuclei with Nq, > 5 is defined as AE^ = 
£^ + 3e^Q, + -Eq and it is suggested that the surface tension energy of the last cluster E^ = E^^'^ . 
The Coulomb energy of the single pn-pair in case of odd Z\ AEpn = 3e^„^ + E^,/;^ and the surface 
tension energy of the pn-pair E^^ = E^^. Therefore, the empirical values of the surface tension 
energy are obtained from 

E'^ = E'J; Ef = E'^ + E'pl^, (6) 

5 

where Ef^ = AE^ — — 3e^„ and E^^^ = AEpn — '^^pna- The values are presented in Table 1. This 
way to obtain empirical values of E'-' and E^^ makes the total binding energy equal to (l)/(2), 
see Table 1. 

The formula to calculate the nuclear binding energy is written as follows 

^ ^nuc ^ + Ean, (7) 

For the hghtest nuclei with Z < 8, < 4, E'*=0. The values of E^'^" and E'^ in this case are 
calculated in accordance with Fig. 1 by adding the corresponding energy portions of the ct-clusters 
and the pn-pair. In case of the nuclei not having a core the binding energy of excess neutrons -Baat 
stays out of consideration. 

For the other nuclei with Z > 9 E^'^'^ is calculated by (5), E*-^ is calculated by a well known 
formula for a spherical charged body with radius Rc, for odd Z nuclei Rci {Rc and Rci are 
determined by (22) on the number of a-clusters in the nucleus, see Section 3) 

3 Z^e^ 

E^^ is calculated by (6) with E^ = JiRq, so Eq.(6) is rewritten as the following 

Est ^ + '£^^Rl.Ef = E^' + jiRlJ2, (9) 

6 

where E^' = 1.026 MeV (see Table 1 for Z = 10, A^^ = 5), 71 = 0.471 MeV/fm^, Rp^ stands for the 
radius of the pn-pair (21), see Section 3. The binding energy of excess neutrons Ean is calculated 
by (13), see Subsection 2.2. 

The values of E^ (7) for the nuclei with N = Z, Ean = 0, are presented in Table 1 in the 
column of 'AE^ry. 

The well known Weizsacker formula [20] is 

The values AE^ (7) depends on and AA^, AE'^^ (10) depends on A and Z. The least 
differences between the values of the Coulomb energy estimated in these approaches are for the 
nuclei with N = Z, see Table 1 [17]. With Z growing the difference increases. Calculation of 
E'> (7) for l^Gd, Eexp = 1338 MeV, gives |^ = 627 MeV and the nuclear force energy of all 
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inter-nucleon interactions E^"^" + E^* + E^n = 1964 MeV, so AE'' = 1 MeV. Calculation of E^^ 
(10) gives ^ = 565 MeV and aA - /3^2/3 + ^ - = 1902 MeV, so AE^^ = 1 MeV. In 

despite of the fact that AE^^ = AE^, these two approaches have different values of the Coulomb 
energy and the nuclear force energy of all inter-nucleon interactions. 

The values of E^ (7) and E^^ (10) calculated for the nuclei with different number of excess 
neutrons are presented in Fig 2. Every nucleus is presented in the graphs and 'E^^^ by two 
dots, one for Eg^p and the other for the calculated value. 
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Figure 2: The values of (7) and the Weizsacker formula i?^^^(10) calculated for the nuclei A(Z) 
indicated in the corresponding graph. 

The widths of the graphs contain both the width of E^-^p distribution on A and the deviations 
between experimental values and calculated ones. For the nuclei with Z > 100 the values calculated 
by the Weizsacker formula have been used instead of unknown experimental values. The width of 
graph 'E''' in the part corresponding to the nuclei with Z > 100 consists of the both the width of 
and the deviations E^^ - E'', which is within a few MeV. The width of the part of 'E^' is a 
little bit bigger than that of the corresponding part of graph where AE^^ = 0. The equal 

widths of the two graphs for the nuclei with Z < 100 clearly show that the deviations AE^ and 
AE'^^ arc close. 

To test the validity of formulas (5), (8) and (9) independently of the energy of excess neutrons, 
the alpha-particle separation energy i?*^^ = E^{A, Z) — E''{A — 4,Z — 2) as well as the deutron 
separation energy E^^^ = E^{A, Z) — E''{A — 2,Z — 1) for the nuclei with even and odd number of 
excess neutrons AA^ has been calculated. In Fig. 3 the calculated values are given in comparison 
with the experimental data. 

The formula to calculate E^^p does not contain Ean neither for even nor for odd excess neu- 
trons. Neither does the formula to calculate E^^^ in case of even number of excess neutrons. The 
deutron separation energy in the nuclei with even number of excess neutrons splits into two lines 
in dependence on whether the nucleus is of even Z or odd Z. It is clearly seen in the left graph of 
Fig. 3. 

In case of odd number of excess neutrons the experimental energy of the single excess neutron 
is bigger for odd Z nuclei than that of even Z nuclei. It is taken into account in the formula 
to calculate the energy of the last excess neutron (14), see Subsection 2.2. Then the formula to 
calculate the deutron separation energy must have —0.2Enn for even Z and -|-0.2£^„„ for odd Z 
nuclei where Enn is the binding energy of the nn-pair the single neutron belongs to. This is in 
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Figure 3: Calculated values of a-particle and deutron separation energy (solid lines) in comparison 
with the experimental values of a-particle separation energies (squares) and deutron separation 
energies (crosses) for the nuclei of /3-stability path and near it. Left graphs are for the even number 
of excess neutrons and the right graphs are for the odd number. 

agreement with the experimental data, see the right graph in Fig. 3. 

One can see a periodical structure in the experimental data, which refers to the shell effects 
connected with the spin-orbit correlations of the last nucleons, so called valent nucleons. The cal- 
culated values are not sensitive to the details. They describe the a-particlc and deutron separation 
energies for all nuclei with an average deviation from the experimental values in 2.5 MeV, which is 
a small percentage of the total binding energy. 

After simplifying Eqs. (5) and (8) one gets E^""" = 42.110Ar« - 26.100 MeV, E^f" = El""" + 
16.884 MeV, and E^ = 1.848iVy^ for odd Z E^ = 1.848(iVa + Q.bf/"^. 

For the nuclei with Z <29 the empirical values of E^ and i?** have been obtained with a good 
accuracy. Table 1. It is good to have a simple function to calculate E^^ instead of the sum (9). It 
can be successfully approximated in accordance with the representation of the nucleus as a core 
plus four peripheral clusters like it is used in calculations of radii (24), see Section 3. Then the 
function consists of the surface tension energy of the core clusters E^^^^ and the energy of four 
peripheral clusters placed at a distance Rp (21) from the center of mass 

E^' = Ef^, + 471 Ef = E^' + 7ii?pi/2, (11) 

where E^l^^ = {N^ — 5.28){Na — 5)^/^ — 11.5 MeV. The approximation is good for the nuclei with 
Z > 22. 

A simple function to approximate for Z >30, > 15, is the following 

E'' = (Ar„ + 1.7)(iV« - 4)2/3; Ef = (Ar„ + 2.8)(iV« - (12) 

The values of \AE^\ of E^ (7) with using E'^ (9), (11) and (12) as weU as the values of {AE^'^l 
for the most of the cases are less than 0.5% of the total binding energy. 

2.2 Excess neutron binding energy 

As it was shown [10] the core starts growing with from the nucleus with = 5. It is suggested 
that the value of the binding energy Ean of the number of excess neutrons AA'" placed in the core 
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and bound in M pairs, AN = 2M + lis calculated as follows: 

M 

Ean = ^ + Cn, (13) 

i=l 

where stands for the energy of ith. nn-pair, e„ is the energy of the single excess neutron and 

en = EnnM+if^ ^ O.lEnnM+i ' (14) 

where '+' stands for odd Z nuclei and '-' for even Z ones. The phenomenological formula takes 
into account that fact that a sum of experimental energies of single excess neutrons of the nuclei 
with odd Z and even Z approximately equals the energy EnuM+i of the M + 1th nn-pair which 
they belong to. In case of AN = 1, M = EnuM+i — Enm- 

The experimental values of Enm ™ MeV are fitted by an equation 

Enm = 22.5 - 1.358f/^. (15) 

It is known that the experimental values of separation energy of nn-pairs have deviations within 
up to 5 MeV determined by what nucleus loses the pair, which certainly refers to the shell effects. 
Despite the relatively big deviation of the values, the empirical values of the energies of all excess 
nn-pairs -£^AAr(exp) known only for the nuclei with N = Z from the following equation: 

EANiexp) = EexpiZ, N + AN) - E,^p{Z, TV), (16) 

are restricted within relatively narrow corridor. For example, for the nuclei with 21 < Z < 29 
the value of the separation energy of two excess neutrons varies within £^2(ea;p) = 21 23 MeV, 
for four excess neutrons E^f^^^p^ = 42 45 MeV, Egi^^^p-^ = 61 ^ 63 and -E'8(ea:p) = 77 80 MeV. 
Therefore, the parameters in (15) have been obtained by fitting both the experimental values of 
nn-pair separation energies known [18] for 27 nn-pairs and the values of -EAiV(ea;p)- The values of 
Ean (13) have been used here to calculate the binding energies of the excess neutrons filling the 
core, which exists only in the nuclei with Z > 10, N^ > 5. It is evident that Eq. (13) is not proper 
for the exotic nuclei like, for example ^He or ^^Li. 



3 Radius of the last proton position, Coulomb radius and charge 
radii of nuclei 

To estimate the value of the radius of the last proton position Rp in a nucleus A the value of AEpn 
is taken as the Coulomb energy of the last proton. The center of mass in the approach is supposed 
to be in the center of the Coulomb field. The energy of the last proton in the nucleus with an 
even Z consists of the energy of its interaction with the other proton of the a-cluster and the 
Coulomb energy of its interaction with the other Z — 2 protons of the nucleus 

AEpn = + (Z - 2)e^/Rp, (17) 

and for odd Z 

AEpn = (£?na ' ^a) + (Z - l)e^ / Rj,l, (18) 

where Rpi stands for the radius of the position of the single pn-pair. 

It is valid only for the nuclei with Z > 10, because the last proton Coulomb energy for the 
lightest nuclei must be calculated with a different function. For example, one can use the spherical 
function with the Coulomb radius Rc like the Coulomb potential at the small radii less than 
Rc in the Shrodinger equation for proton bound state wave function in a Distorted Wave Born 
Approximation . 
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Another way to estimate Rp comes from an equation for the Coulomb energy for the last a- 
cluster. The Coulomb energy of the cluster consists of the internal Coulomb energy of the cluster 

Coulomb energy of its interaction with the three nearby clusters Se^^ and the long-range part 
of the Coulomb energy of its interaction with the other Na — 4 clusters of the nucleus 

AEa = e^ + 3e^^ + 2{Z - 8)eVi?jva-4, (19) 

where Rj^a-A. stands for the distance between the mass center of the remote — 4 clusters and 

the cluster under consideration. The value Rj^a-A is approximated as Rn^-a = l-2Rp by fitting 
the empirical values of AEq,. For the nuclei with odd Z in accordance with the same logic of taking 
into account the long-range Coulomb interaction, one obtains 

AEpn = ^£^a + 7)eV(1.2i?p + Rpi - Rp). (20) 

The third way to estimate Rp comes from the idea that the long-range Coulomb energy must 
be compensated by the surface tension energy E^ . The latter is expected to be proportional to 
R^. Therefore the last member in the sum (19) is taken equal to E'J = 71 i?^ = 2(Z - 8)eV(l-2i?p), 
and the value of 71 = 0.471MeV/fm^ is obtained from fitting the empirical values of AEa- As a 
result, one obtains an equation to calculate Rp in dependence on only. 

= 2.168(iVa -4)^/3_ (21) 

The value of Rpi is calculated by the same formula with using + 0.5 instead of A^^. 

The values of Rp/pi (21) have been used in (17)/(18) to estimate the values of AEpn for the 
nuclei with Z > 23,26,27,28. The values are given in Table 1. The accuracy of several tenth 
McV is determined by the deviation between two possible ways to obtain AEpn by (17)/(18) or 

(19)/(20). 

Finally, the fourth way to calculate Rp comes from the requirement that the surface tension 
energy is expected to be proportional to the number of clusters on the surface of the liquid drop 
E^ = 71 -Rp = 72-A^^^^, where 72 = 1.645 McV is obtained by fitting the values of AE^- The value 
is supposed to refer to the Coulomb radius Rc- Therefore the equation is 

Rc = 1.869(iV„)i/^ (22) 

For odd Z the value of Rci is calculated by (22) with Nai = Na + 0.5. 

The values Rp/pi obtained by different ways differ within 0.3 fm for the nuclei with 16 < Z < 29. 
For the nuclei with 9 < Z < 15 the difference is bigger. For the nuclei with Z > 22, > 11, 
Rc{22) < Rp (21). 

It was shown [10] that the formula to calculate the radii of nuclei 

R = Ra{Na)^/\ (23) 

successfully fits the experimental radii of the most abundant isotopes R^xp with three close values of 
Ra- Using the value of Ra equal to i?4jj(; = 1.710 fm [22] fits the radii of the nuclei with few clusters 
Na = 3, 4. There is no core in these nuclei. In the nuclei with 5 < Na < 12, where the number of 
clusters in the core is comparable with the number of peripheral clusters, the value Ra = 1.628 fm. 
In the nuclei with Na > 12 the radii i?g^p are well described by the radius Ra = 1.600 fm, which 
is the radius of the clusters of the core [10]. 

Varying the number of the peripheral clusters from 1 to 4 (or 5) in accordance with the com- 
prehensive shell model is the point of the alpha-cluster shell model to calculate the root mean 
square radius Rghi [14,17]. Then the formula is similar to that usually used in the framework of 
the single-particle potential approaches with using occupation numbers assumed in the shell model 
[23] in calculations of root mean square radii. 

It was also shown [17] that the number of the peripheral clusters can be fixed at 4 for the nuclei 
with Na > 12, then the root mean square radius R Rghi where R is calculated as follows 
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For odd Z the charge radius Ri is calculated by the equation 



(24) 



{Na + 0.5)Rl = NaR^ + 0.5Rl^, (25) 

where the radius of the position of the single pn-pair Rpi (21) is weighted by factor 0.5. For the 
nuclei with 5 < A^q, < 10 is calculated as Rshh see (15) in [17]. 

The radii calculated in the representation of the alpha-cluster shell model Rghi have deviation 
< >V2= 0.034 fm with (24) (for odd Z (25)) for the nuclei with 9 < Z < 118 [17]. 

A good agreement between the values of R (23) and R (24) (for odd Z (25)) and the experimental 
radii of the most abundant isotopes R^^p [9,22,24] is shown in Fig. 4. From the figure one can see 
that for all nuclei with Z > 10, > 5, the Eqs. (23) and (24) give close to the experimental data 
values. The values as well as the calculated radii of the a-decay nuclei stable to /3-decay [21] with 
83 < Z < 117 are presented in the tables in [10]. 




Figure 4: Left part: Experimental radii of stable isotopes Rexp in dependence on A and Z. On the 
right part: Charge radii R (23) for the 5 < Z < 23 and (24)/(25) for the Z > 24 (soUd line) and 

1/3 

R = l.eOOiVa' fm (dashed line) in comparison with the experimental radii of the most abundant 
isotopes R^^J^j (crosses). The charge radius of one nucleon of a-clusters rch = R/{2Z)^/^ is indicated 
by the lines corresponding to R. 

The values of R calculated for two nuclei ^Li, ^Be are not in agreement with the experimental 
data. These nuclei do not have enough number of bonds between the clusters and the single pn-pair 
to provide a sufficient nuclear density. For the nuclei with 6 < Z < 83 the deviation of R (23) 
and (24)/(25) from the experimental values Rgxp is < >^/^= 0.050 fm. The deviation between 
the two lines in Fig. 4 for the nuclei with 24 < Z < 118 is < >V2= 0.028 fm. The values 
of experimental radii of 125 stable isotopes Rexp [9,22,24,25,26] are also given in dependence of A 
and Z. The radii Rexp are given without indication of their errors, which in most of the cases are 
within 0.004 0.060 fm. It visually approves that fact [10] that both distributions of radii on A 
and Z are well described by the functions ~ A^^^ and ~ Z^^^. 

4 Conclusions 

The approach based on the a-cluster model proposes some formulas to calculate the binding energies 
and the charge radii of the nuclei of the /3-stability path and around it. The formulas have been 
derived on the basis of the idea of isospin independence of inter-nucleon interactions in the a-cluster 
structure. 
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The approach imphes that the nucleus is a dense package of alpha-clusters. The inter-cluster 
distances are determined by the charge radii of the clusters. Some amount of excess neutrons fill 
in the gap between the matter bodies of the a-clusters of the core [10]. The energy of these excess 
neutrons is described by a smooth function on the number of the nn-pairs (13). The formula to 
calculate the binding energy proper for the nucleus with five a-clusters (7) turned out to be good 
for the other nuclei up to the most heavy ones. The same thing one can say about the formula 
to calculate the charge radii R = Vchi'^ZY^^ with the average charge radius of one nucleon of the 
a-cluster structure Vch ~ 1-01 fm for the all nuclei with > 5, see Fig 4. 

The formula to calculate the nuclear binding energy is evidently different from the well known 
Weizsacker formula. These two approaches give different estimations of the total Coulomb energy 
and the energy due to all inter- nucleon interactions, but the values of the total binding energies of 
these approaches are close. To calculate the charge radii both the approaches propose successful 
but different formulas, one is ~ A^^'^ and the other ~ Z^^^. 

A few useful phenomenological formulas have been found in the approach. These are the 
formulas to calculate the root mean square charge radius, the Coulomb radius and the radius of 
the last proton's position in dependence on the number of a-clusters. Besides, the empirical values 
of the Coulomb energy and the surface tension energy with a good accuracy have been obtained 
for the nuclei with N = Z. 

The values of the Coulomb repulsion energy and the nuclear force energy of the inter-cluster 
interaction, which are 1.925 MeV and 4.350 MeV, correspondingly, could be a good test in a 
microscopic description of the light nuclei ^^C and ^^O in the framework of the few-body approach. 
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